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Motivated by the recent measurements of the Bs — Bs 
collaborations, we probe new physics efltects in the Bq 



mass difference from the D0 and CDF 
— Bq mixing within the context of the 

(d) 



supersymmetric grand unified model (SUSY GUT). We find that new physics effects in -Bs(d) — B^ 
mixing lead to the correlated information in the branching fractions of the lepton flavor violating 
decays, which may serve as a test of the SUSY GUT. We also discuss the implication of such new 
physics effects on the quark-lepton complementarity in the neutrino mixings. 
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INTRODUCTION 

The recent measurements for the Bs — Bs mass differ- 
ence from the D0 [ij and CDF [I collaborations given 
by 

17ps-i < AMJ'^P < 21ps-i (90% CL, D0) , 
AMI'^P = (17.77 ±0.10 ±0.07) ps-i (CDF), (1) 

have triggered to probe new physics effects in 6 — > s tran- 
sition. Although the experimental results are consistent 
with prediction of the standard model (SM) dominated 
by the t— quark exchange in the B'^ — Bs box diagram, 
they do not fully exclude all the possible new physics ef- 
fects in Ai? — 2 transitions. Since flavor changing b — > s 
transitions are very sensitive to new physics, it is worth- 
while to probe them through the Bg — Bs mixing phe- 
nomena 3,,!^, £J. Moreover, the combined analysis of 
the Bg — Bs mixing and the Bd — Bd mixing may pro- 
vide clearer hint on the existence of new physics in fla- 
vor changing transitions, due to a possible cancellation of 
hadronic uncertainties. 

Probing such a possibility of new physics in flavor 
changing transitions is the main purpose of this work. 
As a concrete example of physics beyond the SM, we will 
consider the supersymmetric grand unified model (SUSY 
GUT) with heavy right-handed neutrinos where the im- 
print of large atmospheric neutrino mixing may appear in 
the squark mass matrices. In this model large Dirac neu- 
trino Yukawa couplings can induce large off-diagonal mix- 
ing in the right-handed down type squark mass matrix, 
and there can exist possible correlation between quark fla- 
vor changing processes and lepton flavor violating (LFV) 
decays, r — > Here we want to test such new 

physics in b — > s(d) processes in the SUSY GUT, to find 
out the correlated new physics effects in LFV decays as a 
function of branching fractions of r — > m(6)7- We exam- 



ine in detail how new physics effects in i?s(d) ^ ^s(rf) mix- 
ing phenomena are correlated with the branching frac- 
tions of the LFV decays. Comparing our analysis with 
future measurements of the branching fractions of LFV 
decays may serve as a test of the supersymmetric grand 
unified model. Recently, the similar idea has been pro- 
posed in Ref. fsll : The authors have studied a correlation 
between Bs mixing and LFV decay r — > fij by assum- 
ing that there is no new physics effects in Bd mixing and 
some arbitrary mixing term in the slepton mixing ma- 
trix. However, in our approach, we consider not only 
such a correlation but also the correlation among the ra- 
tio of new physics contributions to Bd and Bs mixings 
and the corresponding ratio of the branching fractions of 
LFV decay Br{li — > Ij'^)- Thus, our approach is less 
dependent on arbitrary SUSY input parameters and may 
also reduce the hadronic uncertainties in those processes. 
Moreover, interestingly enough, our analysis can lead to 
an implication of new physics effects in the quark-lepton 
complementarity 0] between the solar neutrino mixing 
angle and Cabibbo angle, as will be shown later. 

In general, the i?° — B^ mass difference is defined 
as AM, = 2\Mi2{Bg)\ = 2\{B°^\H^f,B=^\Bl)\, where 
H^^^ is the effective Hamiltonian responsible for the 
AS = 2 transition, and the SM prediction Q is given by 



Mlf{B, 



G 



127r2 



- Mb, r)^ Bb, fl, (1^; Utfc )^ Lo (x* ) , (2) 



where Gf is the Fermi constant, xt — mi/m^, Lq is 
an "Inami-Lim" function is a short-distant QCD 

correction, and fs, and Bb, are non-perturbative param- 
eters from which main theoretical uncertainties arise. 

Due to the hadronic uncertainties in the SM predic- 
tion for Mf^{Bq), various estimates of the SM values of 
AM™ have come out so far consistent. In order to esti- 
mate the SM values of AM^^, in particular, we adopt the 



2 



following two results for the input hadronic parameters, 
Bb^ s ■ The first one is from the most recent (un- 
quenchcd) simulation by JLQCD collaboration , with 
non-relativistic b quark and two flavors of dynamical light 
quarks. The second one is from combined results, de- 
noted by (HP+JL)QCD. Lacking any direct calculation 
of Bb^ with three dynamical flavors, it has been sug- 
gested to combine the results of Jb^ from HPQCD col- 
laboration with that of Bb^ from JLQCD. Then, two 
numerical results for AM^^ are given by [1] 

AM|M = [0.52l°;?y ps-i JLQCD (3) 
= [0.69 ± 0.14] ps-i (HP-KJL)QCD, 

AM™ = [16.1 ± 2.8] ps-i JLQCD (4) 
= [23.4 ± 3.8] ps-i (HP + JL)QCD. 

The experimental result for AAI^^^ [12] is known to be 

AM^^P = (0.507 ± 0.004) ps-^ (5) 

The mixing amplitude including new physics contribu- 
tions can be parameterized in a model independent way 

as 



Mi2{B,) = Mf^{Bg)[l + Rg] 
= Mf^{B,)[l + r,e^' 

From Eq. ([6]), we obtain 



Mlf{B,) 



1 - 2Re 



( MM) 
KMWiB,) 



(6) 



(7) 



Using 
relation 



Aq, we get the following 



cos aq ± \/ COS^ ffq 



I q — I -I i.g I - 

or (A, -1)' < < (A, + l 



A2-1 



(8) 



On the other hand, the complexity of the mixing ampli- 
tude Mi2{Bq) leads to the CP violation. In Eq. ([S]), the 
CP phase may be composed of the SM and new physics 
contributions as 



arg(l -I- r^e*' 



(9) 



Then, the CP phase arisen from new physics contribution 
is expressed by 



sm 



/NP 



^(1 + Tq coscTqY + (r, sincr,)^ 



(10) 



From the relations Eqs. (|8ll0p . we can extract useful in- 
formation on the new physics effects in the Bq—Bq mixing 
by using the experimental results for AAf^'^P and Uq. 



SUSY CONTRIBUTION TO B° - s" MIXING 



As is well known, supersymmetric standard model 
(SSM) is one of the most motivated candidates for new 
physics beyond the standard model. In supersymmet- 
ric theories, the effective Hamiltonian H^^^ receives 
new contributions through the box diagrams mediated 
by gluino, chargino, neutralino, and charged Higgs. It 
turns out that gluino exchanges give the dominant con- 
tributions [l3|. It is widely accepted that the mass- 
insertion approximation in which only terms with off 
diagonal elements of squark and slepton mass matrices 
are considered is profitable to treat the supersymmet- 
ric contributions to flavor changing neutral current pro- 
cesses. In the mass-insertion approximation, the gluino 
contribution to the amplitude of Bs oscillation is given 
in terms of the ratio of the gluino mass rrig to the aver- 
age squark mass rrig, x = rn^/rri^, and the down squark 
mass insertions between second and third generations, 
{^ab)23 = ("^j )23/'7i|, where A and B stand for left 
(L) or right (R) handed mixing. 

A general expression for = Ml2{B s) / Mf^ {B s) has 
been given in Ref. [l3| as follows: 



Rs = 



ai{mq,x)[{6tL)l3 + {SU)l3] + a2{mq,x)[{S'[ji)ls 

(4L)23]+«3(m,^,a:)[(<5^^)23(4L)23] 

aiimq,xWlL)23{SU)23] , (11) 



where the coefficients |ai| ~ 0{1), \a2\ < {a^l < 
|a4| ~ O(IOO). For example, |ai| = 7.2, jaa] = 129.8, 
jasi = 205.7 and |a4| = 803.8 for to, = 300 GeV 
and X ~ 1 [6]. We also obtain similar expression for 



Rd = Ml2iBd)/M'l^{Ba) by replacing (5^,)23 (a, 6 = 
L.R) with (^ff,)i3 ™ El^- fTl|) . Here, we note that the 
LR{RL) contributions {6'lj^)ij are strongly constrained 
by the measurement of 6 — > 37 decay, and thus those 
contributions should be very small. From now on, we 
will ignore the LR{RL) contributions to Rd.s^ and thus 
Rs,d is approximately expressed by 

Rs ^ ai[(<Sij^^-|-(4^)23]-fa4[(<5ij23(4fl)23], 

Rd ^ ai[(<5iJ?3-|-(4^)?3]-)-a4[(,5iJi3(4fl)i3](12) 

We also note that the values of the corresponding coeffi- 
cients Oi in Rd are very similar to those in Rg . 

It was pointed out in [ij] that large angles in the neu- 
trino sector may imply large mixing among right-handed 
down-type quarks if they are grand unified with lepton 
doublets, and the imprint of the large atmospheric neu- 
trino mixing angle may appear in the squark mass ma- 
trices as a large — sn mixing effect through radiative 
corrections due to the large top Yukawa coupling. 

In SSM with the right-handed neutrino singlets ac- 
counting for the data on neutrino oscillation, it is known 
that large neutrino Dirac Yukawa couplings can induce 
large off-diagonal mixings in the right-handed down-type 
squark mass matrix through renormalization group evo- 
lution. Those mixings can be parameterized as (i5^fl)ij = 



3 



(m^7 where m? is the average riffht-handed 

down-type squark mass. In particular, the large bn — sr 
mixing generated in turn feeds into new physics effects 
in B physics, and thus there may be large enhanced Bs 
mixing. As shown in SUSY SU(5), the off-diagonal ele- 
ments in the mass matrix of the down-type squarks can 
be generated through RG running, and is approximately 
given by Q 



KhJ^^- =^ -^{YjY^hiiml + Al) In 



M< 



GUT 



(13) 



where toq and Aq stand for the universal scalar mass and 
the universal A-parameter for the soft SUSY breaking, 
M* for the scale where the universality of the scalar mass 
is imposed and Mqvt denotes the SU(5) breaking scale. 
Then, (m| )32 and (toj )3i contribute to Bs and Bd 

dRR dan' 

mixing, respectively. 

From the experimental results for AMq and their SM 
predictions, one can extract the allowed regions of the 
parameters \Rq\ and cTq by using Eq. ([8]). Using these 
constraints for \Rq\, we can extract useful information on 
the relevant LFV radiative decays in SUSY GUT context. 



CORRELATION BETWEEN B° MIXING 
AND LEPTON FLAVOR VIOLATION 



Let us discuss how the new physics effects extracted 
from measurements of the Bq mixing can be related to 
the lepton flavor violation in the context of SUSY GUT. 
In the case that the RR contribution to Rq dominates 
over the others in Eq. (jlip . the new physics contributions 
to Rq are approximately given by 



Rs 
Rd 



2 

23' 
2 

13- 



(14) 



Since the term {5'^j^)ij is proportional to {YjY^)ij in 
SUSY GUT, we can obtain the following simple relation. 



Rd 



(YjY,) 



(15) 



13 



Therefore, the origin of new physics effects in the Bq mix- 
ing is from the SUSY seesaw in this case. 

On the other hand, the SUSY seesaw model we con- 
sider can lead to sizable effects on the LFV processes 
such as li — > Ijj due to the new source of lepton flavor 
violation arisen from the misalignment of lepton and slep- 
ton mass matrices, and the branching ratios of the LFV 
decays depend on the specific structure of the neutrino 
Dirac Yukawa matrix. In the context of SUSY GUT, this 
mixing in the charged lepton sector is dictated as same 
as that of the down type quark sector in Eq. ^3]) . As dis- 
cussed in the LFV processes in SUSY GUT models 
can provide a probe of quark-lepton unification. Thus, 
combining the idea proposed in [17| with the analysis 
based on Bq mixing, we can further probe quark-lepton 



unification. The contribution to the branching fractions 
of the LFV decays due to the slepton mass term is roughly 
given by [l5l | 

3 4 

Br{h - In) ^ ^^l(4«)..f tan2/3, (16) 

where ms is a supersymmetric leptonic scalar mass scale 
and we used a rough GUT relation, {5'^R)ij — {^^ll)^- 
We remark that such a GUT relationship for the param- 
eter (5's must be corrected down to the typical mass scale 
of the right-handed Majorana neutrinos Mr because the 
slepton mass term gets additional corrections due to RG 
evolution. However, in general, such RG effects as well 
as corrections from RG evolution down to My/ do not 
significantly modify the GUT relations for off diagonal 
elements {^''£^RR^)ij presented at Mqut^ and, further- 
more, the logarithmic scale dependence is suppressed in 
the ratio of branching fractions of LFV processes. Relat- 
ing Eq. (|15p to the expression for Br{li — > ^^7), we can 
derive the following simple relation 



El 

Rd 



Br{T > fij) 

Br{T > e"f) 



(17) 



Therefore, using the experimentally allowed regions of 
\Rs/Rd\, we can predict the ratio of the corresponding 
LFV processes. 

However, recent work reported that the SUSY models 
with the dominant RR mixing would be disfavored by 
the AMs constraints [6j]. In fact, the LL squark mixing 
receives renormalization group (RG) effects through the 
CKM matrix. The evolution from to the weak scale 
AIw leads to the LL mixings such as 



{f^LLm - -^-2^t ^ts — 2 In——, (18 j 



1 



8^2-^ "' Mw' 



(5^Ji3 - -^Y, Vu In—, (19) 



where Yt is the top quark Yukawa coupling, mo is the 
typical soft SUSY scale. It is known that the RG evo- 
lution from the GUT scale in supergravity scenario in- 
duces {51^)23 - 0.04 - A2 and (<5^i)i3 ~ A^. Then, the 
contribution of the double insertion [{5'li^)23{5jij^)23\ in 
Eq. (fTTjl should not be ignored. In fact, as long as the 
size of {Srr)23 is not greater than 0{X), it turns out that 
the term, a4[{S'[ 23 {Srr) 23)], in Eq. pT|) can dominate 
over the term, ai [(<5^fl;)23)]'^, due to \a4\ ~ 100|ai|. Keep- 
ing only the leading term, ai[{S'l^)i3{S'}^j^)^3], in Rs,d, we 
obtain the following simple relation between the ratio of 
Rs/Rd and the ratio of the branching fractions for LFV 



^ But, there exist highly model dependent cases with large mixings 
in neutrino Dirac Yukawa matrices which may destroy the GUT 
relation for 5's due to large change from RG evolution Those 
cases are not relevant to our work. 
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decays, 



Rd 



Vtd 



Br{T — > /i7) 
Br{T — > 67) 



(20) 



In the following section, we will perform numerical study 
in detail. Note that although we keep both contributions 
proportional to the coefficients ai and 04 in our numerical 
analysis, the relation Eq. (PO)) holds reasonably well. 



NUMERICAL ESTIMATES AND DISCUSSIONS 

Let us begin by examining how to extract the allowed 
regions of and aq based on the SM predictions for 
AM^^. As demonstrated in the introduction, we can es- 
timate the SM predictions for AM^^ by using the results 
of the hadronic parameters from JLQCD collaboration 
and the (HP+JL)QCD collaboration. In this paper, we 
only take the estimation by the (HP-I-JL)QCD collabo- 
ration as given by Eqs. ([3]) and ([5]), from which we could 
probe new physics effects in Bd and B^ mixing more con- 
cretely. As described in Ref. by using Eqs. (6-8) 
and experimental results for AM^'^p, we can obtain new 
physics contributions to Mi2{Bq) which is parameter- 



ized in terms of and a^. 



First, the values of Aq are 



extracted to be 



Arf = 0.75 ±0.30 
A, = 0.74 ±0.18. 



(21) 
(22) 



Using these results for Aq and the relation in Eq. ([5]) , one 
can obtain the constrained regions of the parameters rd^s 
and (Td.s, which are presented in Fig. 1. In addition, the 
parameters rd and ad can be further constrained through 
the allowed values of the CP phase of new physics in the 
Bd mixing (p^^ which is obtained by the experimental 
value of the Bd mixing phase <j)d{= 'I^T^{= Psm) + (f>d^) 
For Bd system, in particular, we have used the known 
constraint on (j)^^ given by Q 



iNPi 
»d Unci 



= -(10.1 ±4.6)" 



(23) 



which further constrains rd and ad through the relation 
of Eq. pO]) . Note that at present there is no constraint 

By using the allowed regions of rd^s and ad.s pre- 
sented in Fig. 1, we can estimate the value of the ra- 
tio Brir — > iJ,'j)/Br{T — > 67) through the relation 
Eq. Although it appears that the effects of the 

phases cr^.s vanish away by taking absolute values in 
Eq. (PO)) . the phase dependence is still imprinted in the 
absolute value \Rq\{— rq), as shown in Fig. 1. Please 
note that the approximate relation Eq. (^0)) holds quite 
well although it keeps only the contributions propor- 
tional to the coefficient 04 in Eq. p^ . In Fig. 2, 
we display the scatter plot of the result for the ratio 
of branching fractions Br{T — > fij)/Br{T — > ej) vs. 
the CP phase of new physics (j)^^ ■ Note that the val- 
ues of the CP phase are determined in terms of 
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FIG. 1: The constrained region of the parameters rd.s and ad,s 
obtained by the method given in Ref. Q and based on the 
SM prediction for AM^^ from (HP-FJL)QCD collaboration. 



rs and as through Eq. (|T0|) . As shown in Fig. 2, the 
value of Br{T — > fij)/Br{T — > ej) decreases as (f^^ 
approaches zero. We note that the current limits on 
LEV radiative decays are B{t ^ 67) < 1.1 x 10~^ and 
B{t ^7) < 6.8 X IQ-s from BaBar Collaboration [H], 
and B{t ^ 67) < 1.2 x 10"^ and B{t nj) < 4.5 x 10"** 
from Belle preliminary report [l9| . Therefore, if both 
LEV radiative decays are observed in near future, we can 
narrowly constrain the CP phase </>^^, and SUSY GUT, 
which we now consider, can be ruled out in case that 
Br[T — > fi'y)/Br{T — > 67) is determined to be below 
10^'^. Inversely, if we determine the size of (j)^^ as well 
as that of r^^d narrowly enough, we can predict the value 
of Br{T — > fi'j)/Br(T — > 67), and thus our approach 
can serve as a test of the SUSY GUT. 

We also calculate the branching fraction for the LEV 
decay, r ^j. In Fig. 3, we present the prediction 
of Br{T — > fij) vs. the CP phase of new physics 
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FIG. 2: The ratio Br(r — > ^i'y)/Br{T — > 67) the CP 
phase of new physics (f)^^ ■ 
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FIG. 3: Br{T — > fij) vs. 0^^. Here we adopt a univer- 
sal scalar mass mo — ms ~ rriq — 300 GeV, ^0 = and 
tan/3 — 10. The horizontal lines are the current upper bounds 
of BaBar [l^ and Belle [3 collaborations, respectively. 



Here we adopt a universal scalar mass, mo — ms — mq — 
300 GeV, Ao = and tan /3 = 10. For the scale A/, in 
Eqs. (HH]) and (HH), we take M, = 10^^ GeV. From Fig. 3 
with the upper bound of Belle, one can constrain —120° < 
< 120°. However, we note that the prediction for 
Br{T — > fij) strongly depends on the SUSY parameters, 
whereas the ratio of the branching fractions is almost 
independent of the SUSY input parameters. 

It turns out that our numerical analysis based on 
Eq. US]) leads to 



i?r(r — > fij) 1 



< 



Br{T — > 67) A 



(24) 



where A(~ 0.22) is a parameter of Wolfenstein 
parametrization of the CKM matrix. We note that this 
upper bound is almost the same as that obtained from 
the case with only RR contribution. The existence of 



the upper bound on 



Br{r *M7) 



may lead to an implica- 



:fr(r ^£7) 

tion on the parametrization of neutrino mixing matrix 
with regard to quark- lepton unification 17[ • As shown in 
Ref. [17| , the so-called quark- lepton complementarity re- 
lation, 0soi+Oc = 7r/4 between the solar mixing angle 0soi 
and Cabibbo angle Oc, can be interpreted as a support 
for the quark-lepton unification, which in turn makes it 
possible to decompose the neutrino mixing matrix into 
a CKM-like matrix and a bi-maximal mixing matrix. It 
has also been shown that a particular parametrization 
of J7pMNS with regard to quark-lepton unification could 
be singled out by examining the ratios of the branching 



fractions Br(li 



Ij^)- It turns out that the result 



given by Eq. may indicate that the neutrino mix- 
ing matrix C/pMNS parameterized by C/t(A)C/bimax with 
CKM-like matrix U (A) and the bi-maximal mixing ma- 
trix J/bimax is preferred. 

In summary, motivated by the recent measurements for 
the Bs — Bg mass difference from the D0 [l| and CDF 
collaborations, we have probed new physics effects in 
the Bq — Bq mixing in the context of the supersymmetric 
grand unified model. We have found that new physics 
effects in i3s(ci) — -Bs(rf) mixing lead to the correlated in- 
formation on the branching fractions of the lepton flavor 
violating decays, which may serve as a test of the SUSY 
GUT. We have also discussed the implication of such new 
physics effects on the quark-lepton complementarity in 
the neutrino mixings. 
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